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I. INTRODUCTION 

Double-/? decay is a process in which a nucleus (A, Z) 
decays to a nucleus (A, Z ± 2) by emitting two electrons 
or positrons and, usually, other light particles 



(A, Z) -> (A, Z ± 2) + 2e T + anything. 



(1) 



Double-/? decay can be classified in various modes ac- 
cording to the various types of particles emitted in the 
decay. For processes allowed by the standard model, i.e. 
the two neutrino modes: 2v(3(3, 2vj3EC, 2vECEC, the 
half-life can be, to a good approximation, factorized in 
the form 



'1/2 



— G^ V \M% V \ ^ 



(2) 



where Civ is a phase space factor and M^ v the nuclear 
matrix element. For processes not allowed by the stan- 
dard model, i.e. the neutrinoless modes: Of/3/3, OvfiEC, 
QvECEC , the half-life can be factorized as 



r l/2 



G 0v \M 0v \ 2 \f{ mil U ei )Y 



(3) 



where Gq v is a phase space factor, Mq v the nuclear ma- 
trix element and /(?7ij, C7 e i) contains physics beyond the 
standard model through the masses rrii and mixing ma- 
trix elements U e i of neutrino species. For both processes, 
two crucial ingredients are the phase space factors (PSF) 
and the nuclear matrix elements (NME). Recently, we 
have initiated a program for the evaluation of both quan- 
tities and presented results for /3~/3~ decay [H-Q- This 
is the most promising mode for the possible detection of 
neutrinoless double-/3 decay and thus of a measurement 
of the absolute neutrino mass scale. However, in very re- 
cent years, interest in the double positron decay, /3 + /3 + , 
positron emitting electron capture, EC/3 + , and double 
electron capture, ECEC, has been renewed. This is due 
to the fact that positron emitting processes have inter- 
esting signatures that could be detected experimentally 
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Q. With this article we initiate a systematic study of 
/3 + /3 + , ECf3 + , and ECEC processes. In particular we 
present here a calculation of phase space factors (PSF) . 
A calculation of nuclear matrix elements (NME), which 
are common to all three modes, will be presented in a 
forthcoming publication [6]. 

Estimates of the transitions rates for /3 + /3 + , EC/3 + , 
and ECEC processes were given by Primakoff and Rosen 
already in the 50's and 60's [7|, |8|. Haxton and Stephen- 
son Q calculated half-lives for /3 + /3 + including relativis- 
tic corrections approximately and some non-relativistic 
calculations were done in the 1980's [M E3- In the 
90 's, this subject was revisited by Doi and Kotani (l2T4T3 | 
who also presented a detailed theoretical formulation and 
tabulated results for selected cases. At the same time, 
Boehm and Vogel [l5[ gave more comprehensive results, 
but without a detailed theoretical description. In these 
papers, results for the PSFs were obtained by approxi- 
mating the positron wave functions at the nucleus and 
without inclusion of electron screening. In this article, 
we take advantage of some recent developments in the 
numerical evaluation of Dirac wave functions and in the 
solution of the Thomas- Fermi equation to calculate more 
accurate phase space factors for double-/3 + decay, EC(3 + 
decay, and double- EC in all nuclei of interest. While in 
the case of /3~/3~ our results (and corrections) were of 
particular interest in heavy nuclei, aZ large, where rela- 
tivistic and screening corrections play a major role, in the 
case of /3 + /3 + our results are of interest in all nuclei, since 
in this case there is a balance between Coulomb repul- 
sion in the final state which favors light nuclei, aZ small, 
relativistic corrections, which are large for heavy nuclei, 
aZ large, and screening corrections, which are large in 
light nuclei due to the opposite sign of (3 + /3 + relative to 
/3 _ /3~. Studies similar to ours were done for single-/3 + 
decay and EC in the 1970's [S Ell- 
in this article we specifically consider the following five 
processes: 

i) Two neutrino double-positron decay, 2f/3 + /3 + : 



(A, Z) ->(A,Z-2) + 2e+ + 2v 



(4) 



ii) Positron emitting two neutrino electron capture, 
2vECP+: 



(A,Z) 



(A,Z 



2v 



(5) 
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iii) Two neutrino double electron capture, 2vECEC: 

(A,Z)+2e- -> {A,Z-2) + 2v (6) 

iv) Neutrinoless double-positron decay, 0v(3 + (3 + : 

{A, Z) (A,Z -2) + 2e+ (7) 



v) Positron emitting neutrinoless electron capture, 
QvECP+: 



(A,Z) 



{A,Z-2) + e^ 



The neutrinoless double electron capture process 
OvECEC cannot occur to the order of approximation 
we are considering, since it must be accompanied by the 
emission of one or two particles in order to conserve en- 
ergy, momentum and angular momentum. It will not be 
considered here. 



II. WAVE FUNCTIONS 



B. Electron bound wave functions 

For electron capture (EC) we use positive energy Dirac 
central field bound state wave functions, 



(11) 



where n' denotes the radial quantum number and the 
quantum number k is related to the total angular momen- 
(8) turn, j K = \k\ — 1/2. For i^T-shell electrons n' = 0,k= — 1, 
IS1/2, while for L/-shell electrons n' — 1, k = —1, 2S , 1 / 2 . 
We do not consider here Lu and L///-shells because 
these are suppressed by the non-zero orbital angular mo- 



mentum, 2P 1 /2, 2P 3 / 2 - The bound state wave functions 
are normalized in the usual way 



V , r l 'K M (i") t V'„' KM (r)dr = 



( r ) + /«',« (r) 



(12) 



dr = 1. 



The key ingredients for the evaluation of phase space 
factors in single- and double-/? decay are the scattering 
wave functions and for EC the bound state wave func- 
tions. The general theory of relativistic electrons and 
positrons can be found e.g., in the book of Rose [Hj]. The 
electron scattering wave functions of interest in /3~ /3~ 
were given in Eq. (8) of [2J. In this article, we need 
the positron scattering wave functions, and the electron 
bound state wave functions. 



A. Positron scattering wave functions 

We use, for (i + decay, negative energy Dirac central 
field scattering state wave functions, 



Vw( r ) 



)~K—K 



(9) 



where XZ are spherical spinors and g K (e, r) and f K (e, r) 
are radial functions, with energy e, depending on the rel- 
ativistic quantum number k defined by k = (l—j)(2j+l). 
Given an atomic potential V(r) the functions g K (e, r) and 
f K (e, r) satisfy the radial Dirac equations: 



dg K (e,r) k e - V + m e c 2 

= --9n(e,r) + Ue,r), 



df K {e,r) 
dr 



e-V- m e c k 

r g K {e, r) + -f K {e, r). 

en r 



(10) 



The potential V appropriate for this case is obtained from 
that for electrons by changing the sign of V (Z into —Z). 
These scattering positron wave functions are normalized 
as the corresponding scattering electron wave functions, 
Eq. (12) of [2j, except for the change in sign in the Som- 
merfeld parameter 77 = Ze 2 /hv. 



C. Potential 

The radial positron scattering and electron bound wave 
functions are evaluated by means of the subroutine pack- 
age RADIAL [l9[, which implements a robust solution 
method that avoids the accumulation of truncation er- 
rors. This is done by solving the radial equations by 
using a piecewise exact power series expansion of the ra- 
dial functions, which then are summed up to the pre- 
scribed accuracy so that truncation errors can be com- 
pletely avoided. The input in the package is the potential 
V. This potential is primarily the Coulomb potential of 
the daughter nucleus with charge Z4, V(r) — Zd{ahc)/r 
in case of (3 + decay and the Coulomb potential of the 
mother nucleus with charge Z m , V(r) = —Z m (ahc)/r in 
case of electron capture. As in the case of single-/? decay 
and electron capture we include nuclear size corrections 
and screening. 

The nuclear size corrections are taken into account by 
an uniform charge distribution in a sphere of radius R = 
TqA 1 / 3 with tq = 1.2 fm, i.e. 



V(r) 



I Zi(ahc) 

±z 



(ahc) ( 



3-(r/fl)- 
2R 



r>R 
r < R 



(13) 



i = d,m. The introduction of finite nuclear size has also 
the advantage that the singularity at the origin in the 
solution of the Dirac equation is removed. 

The contribution of screening to the phase space 
factors was extensively investigated in single-/? decay 
[13,[2l|. The screening potential is of order Vs oc zf^ 3 a 2 
and thus gives a contribution of order a = 1/137 relative 
to the pure Coulomb potential Vc oc Z^a. We take the 
screening contribution into account by using the Thomas- 



3 



Fermi approximation. The Thomas- Fermi function ip(x), 
solution of the Thomas-Fermi equation 



dx 2 



3/2 



with x = r/b and 



b =2 T 



3vr\ 2/3 h 2 



-1/3 



0.8853a Z, 



-1/3 



(14) 



(15) 



where i — d,m and ao is the Bohr radius, is obtained 
by solving Eq. (fl"4"| for a point charge Z% with boundary 
conditions 



<P(0) = 1, 
tp(oo) 



2 



(16) 



for /3 + /3 + decay, 



V(0) = 1, 



ip(oo) 



<p{oo) 



1 



(EC) 
(P + ) 



(17) 



for ECj3 + decay (EC, /? + , respectively), and 

<p(0) = h 
(p(oo) = 0, 



(18) 



for ECEC decay. This takes into account the fact that 
the final atom is a negative ion with charge —2, —1 or 
a neutral ion depending on the mode (f3 + (3 + , ECfi + , 
ECEC, respectively). With the introduction of this 
function, the potential V(r) including screening becomes 



V(r) = ip(r) x 



± 
±Z. 



Zj(ahc) 



(ahc) ^ 



3-(r/R)< 
2R 



r>R 
r <R 



(19) 



i = d, m. This can be rewritten in terms of an effective 
charge Z c g = Zi<p(r) where Z e g now depends on r. In 
order to solve Eq. (fl"4l , we use the Majorana method de- 
scribed in 22] which is valid both for a neutral atom and 
negative/positive ion. The Majorana method requires 
only one quadrature and is amenable to a simple solu- 
tion, the accuracy of which depends on the number of 
terms kept in the series expansion of the auxiliary func- 
tion u(t) of Ref. 22]. The solution is smooth for all three 
boundary conditions. It is particularly useful here, since 
we want to evaluate screening corrections in several nu- 
clei. As an example for the resulting ip(x) functions with 
the boundary conditions presented in Eqs. p^|) -([T5]l we 
show in Fig. [T] results for 78 Kr decay. 
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FIG. 1. (Color online) The Thomas-Fermi functions with 
the boundary conditions of Eqs. (|16[) - (|18[) for 7S Kr decay. 
The dotdashed (gray) curve corresponds to the solution for 
P + fi + decay, the dashed (red) curve corresponds to the solu- 
tion for EC, the dashed (blue) curve corresponds to the solu- 
tion for /3 + decay, and the solid (black) curve corresponds to 
the ECEC. 



D. Solutions 

In order to illustrate the effect of finite size and screen- 
ing we show in Fig. [2] the positron scattering wave func- 
tion for e = l.OMeV, and in Fig. [3] the electron bound 
wave function for the l^i/o and 2 < S 1 / 2 states. Comparing 
Fig. [3]with Fig. 2 of Ref. Fig. 2, one can see that the 
effect of screening is larger than in f3~ /3~ and of opposite 
sign, since the electron cloud decreases the magnitude of 
the repulsive potential seen by the outgoing positrons. 



III. PHASE SPACE FACTORS IN DOUBLE-/? 
DECAY 

In order to calculate PSFs for /3 + /3+, EC/3+, and 
ECEC, we use the formulation of Doi and Kotani 



A. Decays where two neutrinos are emitted 

The 2j//3/3 decay is a second order process in the ef- 
fective weak interaction. It can be calculated in a way 
analogous to single-/? decay. Neglecting the neutrino 
mass, considering only S-wave states and noting that 
with four leptons in the final state we can have angular 
momentum 0, 1 and, 2, we see that both + — > + and 
+ — > 2 + decays can occur. We denote by Q t , where i = 
P+f3+,ECP + ,ECEC, the Q-values of the decay. These 
can be obtained from the mass difference between neu- 
tral mother and daughter atoms, M(A, Z) — M(A, Z — 2) 
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FIG. 2. (Color online) Positron radial wave functions /-i(e, r), (;_i(e,r) (panels a and c, respectively) and <?i(e,r), fi(e,r) 
(panels b and d, respectively) for Zd = 34, e = 1.0 MeV and R = 5.13 fm (vertical line). The notations WF1 (dotted lines), 
WF2 (dashed lines), and WF3 (solid lines) correspond to leading finite size Coulomb, exact finite size Coulomb and exact finite 
size Coulomb with electron screening, respectively. 



as 

Q p+p+ = M(A, Z) - M(A, Z - 2) - 4m e c 2 , 
Qec/3+ = M(A, Z) - M(A, Z - 2) - 2m e c 2 , (20) 
Qecec = M(A, Z) - M(A, Z - 2). 

For the total available kinetic energy one also needs to 
take into account the binding energy of the captured 
electron and thus the total available kinetic energies for 
P + f3+, EC (3+ and EC EC modes are 

T p+fi+ = M(A, Z) - M(A, Z-2)- Am e c 2 , 

T ECf}+ = M(A, Z) - M (A, Z-2)- 2m e c 2 - e b (21) 

Tec ec = M(A, Z) - M{A, Z - 2) - e bl - e b2 . 

The values of M {A, Z)-M (A, Z-2) arc shown in Table|U 
Another quantity of interest in the evaluation of the PSFs 
is the excitation energy En of the intermediate nucleus 
with respect to the average of the initial and final ground 



states. 

A= 1 -W 0+ E N -E I = 1 -[M(A,Z) ^ 
- M (A, Z-2)- 2m e c 2 } +E N -Ei, 

illustrated in Fig. 01 As discussed in Ref. 0], the re- 
sults for PSF depend weakly on the values of the en- 
ergies En in the intermediate odd-odd nucleus, as re- 
marked years ago by Tomoda [23] and as shown explicitly 
in our Ref. 0, Fig. 4. We therefore perform all calcula- 
tions in this paper by replacing En with an average value 
(En) and A = 1.12A 1 / 2 MeV as suggested by Haxton and 
Stephenson 0. The error introduced by this approxima- 
tion is discussed in the following Sect.|lVl We emphasize, 
however, that our calculation has been set up in such a 
way as to allow a state by state evaluation, if needed. 
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FIG. 3. (Color online) Electron bound state wave functions <?o,-i( r )i /o,-i( r ) (panels a and c, respectively) and g% -i(r), 
/? _i(r) (panels b and d, respectively) for Z m = 36 and R = 5.13 fm (vertical line) scaled dimensionless with a factor of 
[47r(m e c 2 ) 3 ]~ 1/2 (ftc/ao) 3/2 ao. The notations WF1 (dotted lines), WF2 (dashed lines), and WF3 (solid lines) correspond to 
leading finite size Coulomb, exact finite size Coulomb and exact finite size Coulomb with electron screening, respectively. 
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FIG. 4. Notation used in this article. The example is for 
106 Cd decay. 



1. 2vf3+l3 + decay 



The formulas for 2i//3 + /3 + decay are exactly the same 
as for 2vj3~ (3~ decay described in [2| where now t\ is the 
energy of the first positron, e\ — e Pl , and ti is the energy 
of the second positrorK €2 = e P2 . We use here the same 
approximations as in [2j, that is to evaluate the positron 
wave functions at the nuclear radius 

5-i (e) = 9-i(e,R) r9 „s 
h(e) = h(e,R) ' (26) 

and to replace the excitation energy En in the interme- 
diate odd-odd nucleus by a suitably chosen energy (En), 
giving 



A=-W a + (E N )-E I . 



(24) 
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Nucleus 



M(A, Z) - M(A, Z - 2)(MeV)* 



P+/3+, EC (3+ and ECEC allowed 



78 Kr 


2.8463(7) 


96 Ru 


2.71451(13)^ 


106 Cd 


2.77539(10^ 


i24 Xe 


2.8654(22) 


^Ba 


2.619(3) 




2.37853(27)^, 




EC (3+ and ECEC allowed 


50 Cr 


1.1688(9) 


58 Ni 


1.9263(3) 


64 Zn 


1.0948(7) 


74 Se 


1.209169(49)^ 


84 Sr 


1.7900(13) 


92 Mo 


1.651(4) 


102p d 


1.1727(36)^ 


112 Sn 


1.91982(16)f 


i20 Xe 


1.71481(125}^ 


144 Sm 


1.78259(87^ 


156 Dy 


2.012(6) 


162 Er 


1.8440(30^ 


168 Yb 


1.40927(25)^ 


174 ffi 


1.0988(23) 


184 0s 


1.453(5s£ 


190p t 


1.384(6) 



ECEC allowed 



36 Ar 


0.43259(19) 


40 Ca 


0.193510(20) 


54 Fe 


0.6798(4) 


i08 Cd 


0.27204(55^ 


126 Xe 


0.920(4) 


132 Ba 


0.8440(10) 


138 Ce 


0.698(10) 


152 Gd 


0.05570(18£ 


158 Dy 


0.284(3) " 


164 Er 


0.02507(12^ 


180 W 


0.14320(27^ 


196 Hg 


0.820(3) 


a Ref. 




b Ref. El 




c Ref. [2|| 




d Ref. [11 




e Ref. El 




f Ref. |3 




s Ref. y|| 




h Ref. || 




1 Ref. El 




j Ref. H 




k Ref. || 




1 Ref. 1351 





TABLE I. Mass difference M(A, Z) ■ 
calculation. 



■M(A,Z — 2) used in the 



The phase space factors are then given in terms of quan- 
tities d m 



(K 



1 



(L N ) = 



ei + UJ\ + (E N ) - Ej 
1 

H ; r . 

e 2 + uj 2 + {E N ) - Ej 
1 

e l +LJ 2 + (E N ) -Ei 
1 

e 2 + lji + (E N ) - Ei 



(25) 



These approximations allow a separation of the PSF from 
the nuclear matrix elements and the condition under 
which they are good have been discussed in @. Apart 
from a narrow region around threshold, where the er- 
ror is ~ 1%, the approximations are good throughout. 
For (3 + (3 + decay we have two integrated phase space fac- 
tors G 2 °J and G^J whose explicit expression are given in 
Eqs. (21)-(28) and (34)-(36) of @. Since the calculated 
single- (3 decay matrix elements of the GT operator in 
a particular nuclear model appear to be systematically 
larger than those derived from measured ft values of the 
allowed GT transitions, and this effect is usually taken 
into account by quenching the axial vector coupling con- 
stant gA,eff = Q_gA, it is convenient to separate it from 
the phase space factors G 2 v Also, it is convenient to 
scale the matrix elements with the electron mass, m e c 2 . 
The phase space factors are then in units of yr^ 1 . From 
these we obtain 
(i) The half-life 



2u 
r l/2 



^ 9 A 



m,c 



(ii) The differential decay rate 



ae Pl at 



(26) 



(27) 



in 



where M% v = g\ 



2 M^l I 



(iii) The summed energy spectrum of the two positrons 



dW 2v 



d(e Pl 



A/- 2t ,ln2 dG ' 2 " 



d(e pi +e P2 )' 



(28) 



These three quantities depend only on G 2 J = G^ v ' 
(iv) The angular correlation between the two positrons 



a ei 



dG 2 °J/de pl 



(29) 
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which depend on both G^J and G%~J . 
following subsections 2 and 3, 



Here and in the 



M. 



(2v) 
GT 



A 



GT 



gy_ 

OA 



M 



{2u) 



A, 



whcr 



M, 



M 



(2u) 



(2u) 
GT 

0+ 



= o 



1 ' T1. fl' 







(30) 



and 



The closure energies 



Aqt and Ap could in principle be different, but in this 
article we take Aqt — Ap = A. 

The phase space factors for 2v[3 + (3 + decay are listed 
in Table [TT] column 2, where they are also compared with 
values found from literature [lj, EH (columns 3 and 4) , 
and in Fig. [5j The values in the literature have been 
converted to our notation by removing factors of g\ and 
(m e c 2 ) 2 . The value for 136 Ce should be taken with cau- 
tion because of the very low Q- value. We also have avail- 
able upon request for all 2v(3 + (3 + nuclei in Table [TT1 the 
single positron spectra, the summed energy spectra and 
angular correlations between the two outgoing positrons. 
As examples, we show the cases of 78 Kr — > 78 Se decay, 
Fig. H and of 106 Cd 106 Pd, Fig. The use of a 
screened potential makes a considerable difference com- 
pared to the results obtained when taking into account 
only the finite nuclear size, as shown in Fig. [5] Note the 
difference between the single positron spectra in Figs. [6] 
and [7] for /3 + /3 + and the single electron spectra in Figs. 
6 and 7 of for (3~ (3~ decay. Note also the difference in 
the scale of Fig. [SJ 10 _29 yr _1 , for 2i//3 + /3 + , as compared 
with the scale of Fig. 5 of 0, 10- 21 yr~\ for 2v$- p~ . 



2. 2vECP + decay 

For the calculation of electron capture processes the 
crucial quantity is the probability that an electron is 
found at the nucleus. This can be expressed in terms 
of the dimensionless quantity [l2| 



S 2 , = 1 
n ' K 47r(m e c 2 ) 3 



ficV 



(31) 



{g b n>A R )) + (fn'J R )Y 



where ao is the Bohr radius ao = 0.529 x 10 _8 cm and we 
use for the nuclear radius R = 1.2A 1/,3 £m. For capture 
from the X-shell n' = 0, k = —1, IS1/2 while for capture 
from the L/-shell n' = 1, k = — 1, 2S\/2- 

Denoting by e p the energy of the emitted positron and 
by eh the binding energy of the captured electron, the 
phase space factor can be written as \\2\ 



G 



EC/3-* 



2A 2 (Gcosg) 4 
3 In 2 167r 5 ft 

EC/3+ + e b- e P 



E B U f 

i=0,l Jm 



EC/3 + 



+ ef, + meC 



[(<;- 1 ( ep , J R)) 2 + (/ 1 ( ep ,i?)) 2 ] 

x {{Kn) + (L N ) + (Kn) (Ln)) lo 1 u>2 Pp ce P dw ide p , 

(32) 

where oji and 0^2 ar e the neutrino energies. Now in the 
definition of (Kn) and (Ln) in Eq. (|25]l . e\ — e e = 
— (m e c 2 — eh) is the energy of the captured electron and 
(2 = £p is the energy of emitted positron. Again sepa- 
rating g\ and the electron mass (m e c 2 ) 2 , the PSF are in 
units of yr _1 . From those, we obtain: 
(i) The half-life 



=Q.<S 




100 120 140 

Mass number 



T l/2 



m,c 



2 M {2u) 



(ii) The differential decay rate 



de n 



de r , 



(33) 



(34) 



where A/21/ 



g\ |m e c 2 Af( 



2;/) 



3. 2vECEC decay 



FIG. 5. (Color online) Phase space factors G^ v in units 
(1CT 29 yr" 1 ). The label "DK" refers to the results obtained 
by Doi and Kotani [T3| using approximate electron wave func- 
tions. The figure is in semilogarithmic scale. 



In the case of double electron capture with two neu- 
trinos the energies of the electrons are fixed and the two 
neutrinos carry all the excess energy. The equation for 
PSF then reads [H: 



8 



Gt^CLO^V" 1 ) 




(io- 2 V _1 ) 


G^ EC {\Q- 2 \jr- 1 ) 


Nucleus This work DK BV 


This work 


DK BV 


This work DK BV 



40 Ca 

54 Fe 
i08 Cd 

126 Xe 
132 Ba 
138 Ce 

158 Dy 
180 W 
196 Hg 



P+f3+, ECP+ and EC EC allowed 



78 Kr 


9770 


13600 


16000 


385 


464 


390 


660 


774 


136 


96 Ru 


1040 


1080 


1230 


407 


454 


350 


2400 


2740 


433 


106 Cd 


2000 


1970 


2420 


702 


779 


652 


5410 


6220 


1120 


124 Xe 


4850 


4770 


5410 


1530 


1720 


1408 


17200 


20200 


3500 


130 Ba 


110 


47.9 


59.2 


580 


549 


420 


15000 


16300 


2590 


136 Ce 


0.267 


0.559 


0.795 


190 


253 


192 


12500 


15800 


2420 





ECP+ and ECEC allowed 






1.16xl0~ 6 


1.05xl0~ 6 


Pi /lOO 


58 Ni 


1.11 1.16 


1.00 


15.3 


64 Zn 


3.81xl0~ 9 


3.83xl0~ 9 


1.41 


74 Se 


1.09xl0~ 5 


8.39xl0~ 6 


5.656 


84 Sr 


0.729 


0.616 


93.6 


92 Mo 


0.206 


0.164 


208 


102p d 


1.62xl0~ 6 


7.16xl0~ 7 


46.0 


112 Sn 


4.95 


4.33 


1150 


i20 Xe 


0.730 


0.524 


888 


144 Sm 


2.49 


1.98 


5150 


156 Dy 


25.3 


20.2 


17600 


162 Er 


6.40 6.69 


5.29 


15000 


168 Yb 


0.00979 


0.00763 


4710 


174 Hf 


l.OOxlO" 9 


2.77xl0~ 14 


1580 


184 0s 


0.0299 


0.0156 


12900 


190p t 


0.00588 


0.00235 


12900 



ECEC allowed 



1.25x10" 
0.0469 
0.0207 
46.1 
39.1 
18.4 
0.183 
0.00156 
821 



17.0 



18100 



0.0887 
3.01 
0.281 
1.08 
17.9 
24.0 
9.14 
235 
173 
982 
3100 
2770 
890 
310 
2240 
2290 



TABLE II. Phase s pac e factors G2v obtained using screened exact finite size Coulomb wave functions. For comparison, values 
of Doi and Kotani [l 2f ] and Boehm and Vogel are also shown. They have been extracted from [l^] and pl| by removing 
g\, and converting to yr _1 units. 



riECEC 

( - T 2i/ 



2i 2 (Gcosfl) 
3 In 2 167r 3 ?i 



-(m e c 2 )' 



E 

*,i=o,i 



B 



B 2 



[{K N ) 2 + (L N ) 2 + (K N ) (ijv)) w 2 w|dwi. (35) 



I 

In this case, in the definition of (Kn) and (Ln) in (i) The half-life 
Eq. P5]l. ei = e ei = — (m e c 2 — e^) is the energy of the 
first captured electron and £2 = £e 2 = -(^eC 2 — e& 2 ) is 
the energy of the second captured electron. The values 
obtained are listed in Table [TT] column 8 where they are 
compared with previous calculations (columns 9 and 10), 
and in Fig. [TO] From G E ^ EC we can calculate: 



'1/2 



_ n EC EC' 4 
— U 2v 9 A 



m e c 2 M^ 



(36) 



While in the case of /3 + /3 + and ECf3 + decay all three 
calculations agree within a factor of ~ 1.5, in the case 
of ECEC decay, the calculation reported in the book 
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e p -m e c 2 (MeV) e pi +e P2 -2m e c 2 (MeV) e p -m e c 2 (Me\) 

FIG. 6. (Color online) Single positron spectra (panel a), summed energy spectra (panel b) and angular correlations between 
the two outgoing positrons (panel c) for the 78 Kr — > 78 Se 2v/3 + f3 + -decay. The scale in the left and middle panels should be 
multiplied by A/21/ when comparing with experiment. In panels a and b the upper, solid curve is obtained when taking into 
account finite nuclear size and electron screening, while the lower, dashed curve presents spectra obtained when taking into 
account only the finite nuclear size. In panel c these two calculations coincide. 




e p -m e c 2 (MeV) e Pi +e p -2m e c 2 {MeN) e p -m e c 2 (MeV) 

FIG. 7. Single positron spectra (panel a), summed energy spectra (panel b) and angular correlations between the two outgoing 
positrons (panel c) for the 106 Cd — > 106 Pd 2i//3 + /3 + -decay. The scale in the left and middle panels should be multiplied by A/21/ 
when comparing with experiment. 



of Boehm and Vogel [15( , disagrees with other two by a 
factor of approximately 4. The origin of this discrepancy 
is not clear. The values in Table [TT1 have been converted 
to units yr _1 using the same procedure in all three cases, 
EC/3+, and EC EC. Since apart from the factor 
of 4, the behavior with mass number of G^ EC in 
is the same as in the other two calculations, it may be 
simply due to a different definition of G E ^ EC ■ Note that 
the scale in Fig. M 10~ 24 yr _1 , for IvECEC is very 
different from that for 2i//3+/3+, 1(T 29 yr" 1 , in Fig. |5] 
due to a much larger Q-value 



B. Neutrinoless modes 

As discussed in Ref. [4|, several scenarios of neutrino- 
less double beta decay have been considered, most no- 
tably, light neutrino exchange, heavy neutrino exchange, 
and Majoron emission. After the discovery of neutrino 
oscillations, attention has been focused on the first sce- 
nario and the mass mode, where the transition operator 



is proportional to (m„) /rn e . In this article we present 
phase-space factors for the mass mode. Phase-space fac- 
tors associated with the other modes, called (A) and (rj) 
in Ref. (23|, will form the subject of a subsequent publi- 
cation. 



1. 0u/3 + /3 + decay 

The equations for 0^/3 + /3 + decay are exactly the same 
as for 0v(3~ (3" decay described in |2] where now t\ is the 
energy of the first positron, e\ — e Pl , and £2 is the energy 
of the second positron, £2 = £p 2 ■ There are also here two 
quantities G^J and G^J in units of yr" 1 from which one 
can obtain: 
(i) The half-life 

"t^F 1 = GgXfli ^ 2 M<°"> 2 . (37) 



10 



1000 




100 120 140 

Mass number 

FIG. 8. (Color online) Phase space factors in units 

(1CT 24 yr" 1 ). The label "DK" refers to the results obtained 
by Doi and Kotani using approximate electron wave func- 
tions. The figure is in semilogarithmic scale. 

The obtained PSFs are listed in Table [TT1 column 5 
where they are compared with previous calculations 
(columns 6 and 7), and in Fig. [5J The very small values 
in the second part of the table should be taken with 

caution in view of their very small Q- value. 
An example of single positron spectrum is shown in 
Fig. E This figure is for 106 Cd ^ 106 Pd 2uECP+ -decay. 



8.0x10 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 O.i 

e p -m e c 2 (MeV) 

FIG. 9. Single positron spectra for the 106 Cd -> 106 Pd 
2vEC/3 + -decay. The scale should be multiplied by Miv when 
comparing with experiment. 



(ii) The differential decay rate 



dc 



in 



dc,, 



(38) 



where.^ = 9 a \M^\ 2 . 

Both the half-life and the differential decay rate, are given 
in terms of G^J = Gq„ 13 . 

(iii) The angular correlation between the two positrons 

dG$/de Pl 



a(e Pl ) 



dG { °J/de pi 



(39) 
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FIG. 10. (Color online) Phase space factors G^y EC in units 
(10~ 24 yr" 1 ). The label "DK" refers to the results obtained 
by Doi and Kotani P3] using approximate electron wave func- 
tions. The figure is in semilogarithmic scale. 



G: 



/3 + /3 + (10 -2u 



Nucleus This work DK BV This work DK BV 
P+/3+, EC /3+ and EC EC allowed 



78 Kr 


250 


293 


59.4 


6.37 


7.11 


96 Ru 


84.5 


90.7 


12.2 


9.62 


10.8 


i06 Cd 


96.2 


102 


14.5 


13.0 


14.7 


124 Xe 


114 


123 


18.1 


19.7 


22.9 


130 Ba 


25.7 


21.0 


1.67 


17.6 


19.8 


136 Ce 


2.42 


3.55 


0.175 


15.3 


18.7 



50 Cr 
58 Ni 

64 Zn 
74 Se 
84 Sr 
92 Mo 

102p d 

112 Sn 
i20 Te 

144 Sm 

156 Dy 
162 Er 
168 Yb 
174 Hf 

184 0s 

190p t 



ECP+ and ECEC allowed 
0.0887 

1.21 1.30 
0.0507 
0.230 

1.94 

1.92 
0.287 

5.20 

3.92 

8.11 

15.2 

12.9 15.8 
4.23 
0.0272 
7.04 
5.57 



TABLE III. Phase space factors Gov obtained using screened 
exact finite size Coulomb wave functions. For comparison 
value of Doi and Kotani [13j an Boehm and Vogel [l 51 ] are also 
shown. These are extracted from [p| and [ill by removing 



g A and converting to yr units. 



The values of Gq v are shown in Table IIIII column 2 
where they are compared with previous calculations (col- 



11 



umn 3 and 4), and in Fig.QT] In this case, our calculation 
and that of [131 ] disagree with the calculation reported in 
15;] by a larger factor. 
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Mass number 

FIG. 11. (Color online) Phase space factors Gq u 13 in units 
(10~ 20 yr -1 ). The label "DK" refers to the results obtained 
by Doi and Kotani [12| using approximate electron wave func- 
tions. The figure is in semilogarithmic scale. 

We also have available upon request the single electron 
spectra and angular correlation for all 0^/3 + /3 + nuclei in 
Table HTH An example, 106 Cd decay, is shown in Fig. IT2l 



dius, R. We have taken for the atomic mass M(A, Z) the 
available experimental values with errors shown in Ta- 
ble |U The more accurate values used in this table account 
for some of the differences between our calculated values 
and those of [HI, EH EH obtained with older values of the 
atomic masses. We estimate the error here as a multiple 
of SQ/Q, where SQ is the error in Q. The nuclear ra- 
dius enters the calculation in various ways, most notably 
the evaluation of the positron wave functions at the nu- 
cleus g-i{R) and f x {R) and for ECP+ and ECEC, the 
electron probability, B 2 , K . We have taken R — tqA 1 ^ 3 
with ro = 1.2fm. An estimate of the error here is ob- 
tained as in single-/?" 1 " decay and single EC [2l[ by ad- 
justing ro for each nucleus, A, Z, using the experimental 
value (r 2 }exp from electron scattering. Finally, another 
uncertainty is introduced by the average excitation en- 
ergy, (En)- In [2| we estimated this uncertainty which 
affects only the 2v processes by comparing the results of 
the calculation with (E N ) = 1.121 1 / 2 MeV with that of 
the single state dominance model, for example in 106 Cd 
decay, the ground state of the intermediate 106 Ag nucleus 
is 1 + , giving En = E ± + = O.OMeV. The resulting error 
is of the order of few percent. Screening corrections play 
a minor role and do not introduce a large error. The 
situation is summarized in Table IIVI 



2. 0uEC/3+ decay 

In case of neutrinoless positron emitting electron cap- 
ture, the energy of the emitted positron is fixed and the 
equation for PSF reads }13| : 

G ECfi+ _ 1 2 (Gcosfl) 4 2 2 5 

x E B l-i [(.9-i (e P ,^)) 2 + (fi(t P ,R)) 2 ] P P ce p . 

i=0,l 

(40) 

The PSF are in units yr _1 , and from those we can cal- 
culate: 

(i) The half-life 



T l/2 



_ r EC[3+ 4 
— b 0i/ 9 A 



(m v ) 


2 











(41) 



The values obtained are listed in Table IIIII column 5 
where they are compared with previous calculations (col- 
umn 8) and in Fig. [13] In this case only calculations of 
1131 are available. 





Q-value 


10 x SQ/Q 




Radius 


1.0% 




Screening 


0.10% 




(E N ) 


model dependent (< 1%) 


0^/3+/3+ 


Q-value 


4 x SQ/Q 




Radius 


9% 




Screening 


0.10% 








2vEC/3+ 


Q- value 


8 x SQ/Q 




Radius 


1.0% 




Screening 


0.10% 




{En) 


model dependent (< 1%) 


0isEC/3+ 


Q-value 


2 x SQ/Q 




Radius 


9% 




Screening 


0.10% 




{En) 




2vECEC 


Q- value 


6 x SQ/Q 




Radius 


1.0% 




Screening 


0.10% 




{En) 


model dependent (< 1%) 



IV. ESTIMATE OF THE ERROR 



TABLE IV. The estimate of uncertainties introduced to phase 
space factors due to different input parameters. 



Two main sources of error in the evaluation of the 
phase space factors are the Q- values and the nuclear ra- 
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FIG. 12. Single positron spectra (panel a), and angular correlations between two outgoing positrons (panel b) for the 106 Cd 
— > 106 Pd Of/3 + /3 + -decay. The scale in the left should be multiplied by Nq v when comparing with experiment. 




modes, but in Ref. [37j this value is believed to be for the 
2vECEC process, because other modes are strongly sup- 
pressed. Our calculations in Table [TT1 support this state- 
ment and we thus take T^cec = ^.2 ± 0.5) x 10 21 yr 

in 130 Ba, from which we extract the value of \M^ ECEC \ 
in Table El 

We note that the value we extract is comparable to the 
values extracted from 2vj3~ j3~ decay in Q ■ To emphasize 
this point we show in Fig. [14] the newly extracted value 
in comparison with all others. 



100 120 140 

Mass number 



FIG. 13. (Color online) Phase space factors in units 

(10~ 18 yr~ x ). The label "DK" refers to the results obtained 
by Doi and Kotani pl| using approximate electron wave func- 
tions. The figure is in semilogarithmic scale. 



USE OF PHASE SPACE FACTORS 



The main use of PSFs discussed in this paper is to 
calculate half-lives for the EC (3+ and EC EC de- 

cay, by combining them with a calculation of the NME, 
the only constraints being that the NME are defined in 
a way consistent with Eq. ([26]), 2^/3+/3+, Eq. (|33| . 
2vECf3 + , Eq. (JH, 2vECEC, Eq. ([37]), 0i//3+3+, and 
Eq. gl]), QvEC(3 + . The calculation of half-lives with 
matrix elements obtained from IBM-2 will be presented 
in forthcoming publication [6j. Here we use the cal- 
culation of PSF to extract the dimensionless quantity 
g\\{m e c 2 )M 2 v\ 2 = |A/|5| 2 as done in the case of f}~ f3~ 
decay reported in [2|. 

For two neutrino double positron decay and compet- 
ing modes the only positive experimental half-life re- 
sult is from a geochemical experiment in 130 Ba [36j j : 
lfj 2 = (2.2 ± 0.5) x 10 21 yr- In geochemical ex- 
periments, it is not possible to disentangle the different 
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FIG. 14. (Color online) Effective nuclear matrix elements 
|M|J| extracted from the experimental 2^/3/3 half- lives as a 
function of mass number. The red triangle is the value ex- 
tracted from t^2 CEC ■ F° r comparison the grey dots and 
squares represent the values extracted from f3~ f3~ experi- 
ments. 



VI. CONCLUSIONS 

In this article, we have reported a complete and im- 
proved calculation of phase space factors for 2i/f3 + f3 + , 
2vECfi + and 2vECEC, as well as 0^+/3+ and 
0isEC(3 + double-/? decay modes, including half-lives, 
single positron spectra, summed positron spectra, and 



13 



Nucleus 


Gf l/ CBC (10' 21 yr- 1 ) 


r 2,£CBC (10 21 yr) exp a 


1 M cS 1 


130 Ba 


15.0 


2.2 ±0.5 


0.174 ± 0.017 



a Ref. [H 

TABLE V. Experimental 2vECEC half-lives and the corresponding effective nuclear matrix elements \M^, ECEC 



positron angular correlations, to be used in connection 
with the calculation of nuclear matrix elements. Apart 
from their completeness and consistency of notation, we 
have improved the calculation by using exact Dirac wave 
function with finite nuclear size and electron screening. 
The program for calculation of phase space factors has 
been set up in such a way that additional improvements 
may be included if needed (P-wave contribution, finite 
extent of nuclear surface, etc.) and that it can be used 
in connection with the closure approximation, the single 



state dominance hypothesis and the calculation with sum 
over individual states. 
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